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SUR 

L'ÉVALUATION APPROCHÉE 

DES AIRES PLANES 



Objet du présent mémoire. 



1. Formuks anciennes. 11 existe un grand nombre de formules 
pour le calcul approché des aires planes. La plus ancienne, dont 
l'idée première est due à Côles et que Newton a insérée dans 
les Principes (Livre ïll, lemme V), repose sur la substitution, 
ft la courbe donnée, d'une parabole de degré n, ayant avec 
celte courbe (» + 1) points communs. L'expression approchée 
de l'aire ne contient que les {n -h 1) ordonnées communes à la 
courbe et à la parabole et les distances qui les séparent. Gauss 
a prouvé (Werke , IIÏ, 165-196, 202-206) qu'en choisissant 
convenablement ces distances, on pouvait obtenir une formule 
rigoureusement exacte, dans le cas où la eourbe donnée est une 
seconde parabole de degré (2« — 1) ('). Malheureusement, ni 
Gauss, ni aucun de ses continuateurs n'a cherché, que nous 
sachions, la valeur maxima de l'erreur que comporte sa for- 



i fioatielle correspondance rnalhémalique , t. VI, pp. 396-402, une 
^s simple de ce théorème de Gauss, démonstralien due à M. Ë, Galalan. 
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inule, dans le cas où la courbe dont on veut obtenir Taire appro- 
chée est une ligne continue quelconque, au lieu d'être une 
parabole de degré (2w — 1). 

La formule sommatoire de Maclaurin et d'EuIer n'a pas cet 
inconvénient ; on est parvenu à trouver la valeur de l'erreur 
maxima k laquelle son emploi peut conduire (*). Mais elle pré- 
sente d'autres désavantages. L'aire cherchée, pour une courbe 
y = fx, est exprimée, non seulement au moyen d'un certain 
nombre de valeurs de la fonction f, mais aussi des dérivées suc- 
cessives de cette fonction. En outre, l'erreur qu'elle comporte 
ne peut être calculée d'une manière approchée qu'en cherchant 
entre quelles limites varie une certaine dérivée de la même 
fonction. Par suite, la formule de Maclaurin et d'Euler est, en 
pratique, d"im emploi extrêmement peu commode. 

On peut déduire de la formule de Maclaurin et d'Euler, direc- 
tement ou indirectement, un grand nombre de formules plus 
simples et plus pratiques : la formule des trapèzes, la (première) 
formule de Simpson, la formule de M. Catalan, la formule de 
Poncelet, la formule de M. Parmentier, etc. La (première) for- 
mule de Simpson, qui semble de beaucoup la plus exacte de 
toutes celles que nous venons de citer,, peut s'obtenir aussi 
comme cas particulier de celles de Côtes, de Newton ou de 
Gauss. 

2. Objet de ce mémoire. La démonstration habituelle des for- 
mules de Poncelet et de M. Parmenlier ne laisse rien à désirer 
au point de vue de la rigueur ou de la simplicité et elle conduit 
tout naturellement à une estimation approchée de l'erreur com- 
mise en les employant. Mais II n'en est p"as de même pour les 
autres formules, qui, au fond, sont établies d'une manière à peu 
près empirique. Récemment, en étudiant les écrits de M. le 
général Parmentier sur les quadratures approchées, nous nous 
sommes aperçu qu'on pouvait démontrer toutes les formules 



lir, par exemple, Houel, Catt 
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énumérées plus haut avec la même facilité que celle de Poncelet, 
estimer l'erreur maxima correspondante et la représenter géo- 
métriquement, souvent de plusieurs manières. De plus, on peut 
trouver un grand nombre d'autres formules, dont une au moins 
(dueàM.Ch. Dupain), deux peut-être, semblent avoir une valeur 
pratique, à cause de leur simplicité et de leur exactitude. 

Nous nous proposons d'exposer ici, d'une manière systéma- 
tique, l'ensemble des recherches auxquelles nous avons été 
conduit. La partie la plus élémentaire de celte étude a déjà 
paru dans la Mathesis (t. 1, pp. 17-22, 55-36, février et mars 
1881). 

3. Principaux résultats obtenus. Parmi les résultats obtenus, 
nous signalons spécialement, au point de vue pratique, l'estima- 
tion graphique de l'erreur que comportent les diverses formules 
et particulièrement celle des trapèzes. 

Si l'on représente par AD l'erreur maxima pour cette dernière 
formule, les erreurs analogues, pour les autres formules, sont 
respectivement ; 

i 

-kÙ pour ta formule de Poucelcl; 

-hh poiii'la preiiiicreformulcdcSimpsoii, et pour 
celle de M. Parmenlier et la formule de 
M. Dupain {n° 7); 

-AD pour la seconde formule de Simpson et |ioni' 
la formule de M. Catalan; 

5 

- hU pour les formules de Poncelet, de M. l'armen- 

tier et de Simpson, quand on les étend au 

cas d'un nombre pair d'ordonnées (n" 21); 

4 

-hD pour la formule nouvelle fn' 15), 
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Ces limites de l'erreur peuvent servir, dans un cas donné, à 
voir rapirfement, sur la figure même, par une construction plus 
simple que celle que Poncelet a indiquée pour sa formule, quelle 
grandeur il convient de donner à h pour obtenir une approxima- 
tion déterminée d'avance. 

Une fois h choisi, les calculs qui servent à trouver une valeur 
approchée de S, par l'une quelconque de ces formules {celles de 
Poncelet et de M. Parmentier exceptées), donnent en même temps 
une nouvelle limite de l'erreur commise et cette nouvelle limite 
est inférieure à celle qui vient d'être indiquée, même pour la 
formule des trapèzes. 

Voici encore deux résultats sur lesquels nous attirons l'atten- 
tion du lecteur : 

1" Quand le nombre des ordonnées est pair, la meilleure for- 
mule pratique est celle de M. Catalan , qui est, pour ce cas, la 
transformée de celle de Simpson , et qui peut d'ailleurs être 
employée même si le nombre des ordonnées est impair (n" 21). 

2" La formule de M. Dupain (n" 7), presque aussi simple que 
celle de Poncelet, a une exactitude presque égale à celle de 
Simpson , et est aussi exacte que celle de M, Catalan. 



Préliminaires. 

4. Estimation de l'erreur commise dans (es approximations. 
I. Si une quantité inconnue S est plus grande que m, plus petite 
que M, on peut dire que S = m, avec une erreur, par défaut, 
plus petite que M — m, ou que S = M, avec une erreur, par 
excès, plus petite que M — m. 

II. Si l'on prend S = j (M ■+- m), l'erreur commise est tout 
au plus égale à ^(M — m), mais on ignore quel en est le sens. 

III, Si l'on sait que S est plus rapproché de M que de m, 
S étant alors supérieur à ^ {M + m), on peut prendre pour S 
une valeur intermédiaire entre M et ~ (M ■+■ m), par exemple 
l (2M -+- m). En posant S ^^ | (M -î- m), l'erreur sera par excès 
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si s est compris emre ^ (M + m) et ^ (2M ■+- m), et elle sera 
moindre que 

g(2M + m)-i(M -t- m) = i(M -m]. 

L'erreur sera par défaut si S est compris entre M cl f (2M + m), 
et elle sera moindre que 

i i 

M (2M ^-m) = -(M — m). 

Donc, en tout cas, si S est compris entre ^ (M -i-m) et M, et 
si l'on pose S=^(2M+«i),!'erreur est moindre que \ (M — m). 

IV. En général , si y. est une quantité comprise entre m et M , 
en faisant S=fx, on commet une erreur plus petite que la plus 
grande des deux quantités fi — m, M — p. 

5. Données et notations. Considérons (Fig. 1, 2 ou S) une 
aire S, comprise entre un arc de courbe AL, une droite fixe al 
et les perpendiculaires Aa, Ll abaissées sur celte droite des 
extrémités de l'arc AL. Supposons, pour fixer les idées, que 
l'arc AL soit, dans toute son étendue, concave vers la droite 
al (*). Divisons la base al en un nombre pair de parties égales, 
en dix, par exemple. Par les points de subdivision, b, c, d, e, f, 
g, h, i, k, élevons des perpendiculaires à «/. Désignons ces per- 
pendiculaires ou ordonnées «A, 6B, ... Ih, dont les extrémités 
A, Bj ...L sont sur la courbe, par t/i,y2) •■-' J/iii soit A la distance 
connmune de deux quelconques de ces ordonnées. Nous repré- 
senterons par E, la somme des deux ordonnées extrêmes, par I, 
la somme des autres ordonnées de rang impair, par P, celle des 
ordonnées de rang pair. De celte manière, 

I=!/î-i-.'A+y7 + yo, 

p = y, -l- j,, -t- y. -t- 1/3 -4- y,„. 



ES quantités rf. rf', d",j', se 
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Nous posons ensuile 

-1 a 

2 2 

(("=2^— rf', 
5 = 2P — E — 21. 

■II 

Formules de Pohcelet, de M. Parmentier et de M. DupAiN. 

e. Première limite supérieure et 'première limite inférieure 
de S. I. Par les exirémités B, D, F, H, K des ordonnées de rang 
lair (Fig. 1), menons à !a courbe des tangentes lerminées enB, 




ei Ba, D, et Dï.F, ei F,, H) ei Ilg, K, elK^, aux ordonnées 
voisines. La somme des trapèzes 

uB,B,c, cD.Die, eV.V^g, ïH,H,»', (K,K,(, 



yGoosle 



est supérieure à S. Ces trapèzes ont pour mesure 

2A!/î, 2%4, 2%„, 2%8, 2%i„. 
Donc, 

M = ft {2^3 + 2j/, -H 2)/8 -H 2^8 + 2!/io) = 2/iP. 

il. Le polygone aABDFHKL/ a une aire m inférieure à S. 
On a d'ailleurs 

"1=^ 2(3/'-»-ys)+(y»+y*)+(y»-*-y«)-'-(î'«+y9)+(y»-+';'«)+2(3/<o-i-î/i>) 

= ft[2y, + 2y.H-2y, + 2ya+2;/.„]-Ari{y,+.v,„)-^[y.+y„)l 

ou, d'après les notations du numéro précédent, 

m = 2/(P — Arf, 

JI[. Il est impossible de démontrer d'une manière générale 
que S est plus rapproché de M que de m. Néanmoins, dans le cas 
de courbes à courbure pett. prononcée et lorsque \i est suffisam- 
ment pelil, on reconnaît aisément , en faisant un tracé graphique, 
que l'on a souvent le droit, au point de vue pratique, de sup- 
poser S plus rapproché de M que de m. 

î. Formules de Poncelet, de M. Parmenlier et de M. Dupain; 
interprétation géométrique de l'erreur. I. Puisque S est com- 
pris entre M et m, on peut, d'après le principe II du n" 4, 
prendre approximativement M = p, p étant la demi-somme de 
M et de m; l'erreur s sera inférieure à la demi-différence. On 
aura ainsi la formule de Poncelet : 

S =p (approximativement), 
p = - (M -t- m) = 2/iP — - hd, 

.<i(M-m) o« iftrf. 



yGoosle 



II. Si l'on prend p' = ^ (2M ■+- m), en faisant 

S = p' (approximativcmeni) , 
on aura la formule de M. Parmentier. D'ailleurs, 

5 

L'erreur maxtma est égale à la plus grande des deux quantités 
M — p'=iAd, et p'— M*=§Arf, donc, égale à | hd. 

Très souvent, d'après le principe III du n" 4, et la remarque 
faite au numéro précédent, l'erreur sera inférieure à | (M — m), 
c'est-à-dire à ~hd. 

III. Nous verrons plus loin pourquoi M'. Parmentier a choisi 
pour p' la valeur ^ (2M -t-m), plutôt que i (3M -+■ m), ou toute 
autre valeur comprise entre M et ^ (M + m). Nous verrons en 
même temps qu'une formule préférable à celle de Poncelet et à 
celle de M. Parmentier est la suivante, due a M. Dupain, 

S := TTj (approximativement), 

i hd , 1 , i hd 

T. = »' = 2AP lia — — , 

' 2w - 1 5 3 2n — 1 ô 

2/(rf 1 hd 

' o 2« — ) T ' 

2n étant le nombre des subdivisions h de la base al, donc, dix 
dans le cas de la figure. 

IV. La quantité d et, par suite , les produits hd, ^hd, etc., sont 
susceptibles d'une interprélation géométrique simple. Appelons 
P, Q, les points d'intersection des droitesAL,BK avec l'ordonnée 
moyenne Ff. On aura évidemment 

lis,+s«)=ff, 5 (s. + »..)=«/; 

Par conséquent, 

d = Q/-_Pf_PQ. 
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Il est donc facile de tracer sur la figure même un petit rec- 
tangle PQQ'P' dont Taire est égale à M. 

Par suite, on peut estimer, sur la figure même, avant d'entre- 
prendre aucun calcul, la grandeur des erreurs théoriques |Arf, 
I hd, etc. Si ces erreurs, pour une valeur de h choisie d'avance, 
étaient trop considérables pour le but que l'on se propose en 
cherchant S approximativement, on pourrait recourir à une 
valeur plus petite de h, et ne commencer à estimer numérique- 
ment p, p', ou TTj , qu'avec la certitude de trouver ainsi une valeur 
suliisanimcnl approchée de S (*). 



FomuULlî DES TRAPÈZES ET FORMULE DE SlMPSON. 

Première démonstration. 

s. Deuxième el troisième limite inférieure rfeS. I. Le polygone 
«ACEGIL/ (Fig. 2) a une aire m' inférieure à S ("). On trouve, 
sans peine, 

m'=A[(y, + ys) -H (y, -H y,) + (y, + 2/>) -^(^/7^-y»)-+-i'^-l-^/^^)]- 
0H encore 

m' = A {E -(- 21). 

II. Le polygone oAÎÎCDEFGHIKL/ (**') a aussi une aire t 
inférieure à S. L'aire ( est donnée par la formule 
\ 

' = -/'[(!/( -i-y^} + (.Vî-+- Ih)-^ -^ (.'/a-*-yit>) + (ywH-.'/ii)]> 



-A(E-f-21 +2P). 



r) Ce paragraphe ne coulient absolumeiil rien de neuf el doil élre regardé comnK 
une inlroductioa aux suivants. 

(■■) On ne eoiiçoit pas comment Poiieelet, à qui apparlient la cojisidéralion des aires 
M etin, n'ait pas songea ce polïgone m'quivanaus donner si simplement la démonstra- 
li-ation de la formule de Simpson. 

('") Non tracé sur la figure. 
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», Formule des trapèzes; formule de Simpson. I. La figure 
ou les formules donnent immédiatement 



Donc, d'après le principe II, n» 4, en posant (formule des 
trapèzes) 

S = ( (approxiiiiativemenl), 

on commettra une erreur maxima plus petite que ~ (M —m'). 
Or, 

M — m' = A(2P — E — 21) = hS. 

L'erreur maxima, par défaut, que comporte la formule des 
3S, est donc plus petite que ~ hri. 




II. Puisque S est compris entre M et ( = ^ (M + 7»'), on 
pourra prendre 

S =! s (approximative me ni), 

s étant égale à { (21V! + ?«'). trest la formule de Simpson. On a 
d'ailleurs 

sr=-fe(E + 2I-»-4P), 
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et d'après 4, III, 

j < „ (M — m] ou r 1^^- 

On peut encore écrire s et hS sous la forme 

s=i(M-H2(), /^J=2[M~(). 

On verra plus bas pourquoi la valeur la plus convenable pour s, 
quantité qui doit être choisie entre M et m', ou M et (, est celle 
qui est donnée ici. 

lO. Interprétation géométrique de S. Autres limites de l'erreur. 
I. On peut donner, de 5 aussi, une interprétation assez simple. 
Soient R, S, T les pointa où CI, DH, EG rencontrent F/". On 
aura : 

yi^m^^f, bsL^î^sf, y-i^yj-^Tf. 

Or, ta valeur de <5 peut s'écrire successivement : 

fc(%, + %,+ %, -*- 2^8+ 2j/,o - y, - y„-2y,- 2i/,-2;/,-2y.) 

\^ 2 2 2 2 2 / 



Par conséquent, 

J= 2 [{F/-- T/^ - (T/-- Sn + (S/- R/-) ~ (R/- Q/) + (Q/- P/}] 
= 2[FT-TS+SR-RQ-hQP], 

quantité que l'on peut construire sur la figure. On peut donc 
aussi représenter hd par un petit rectangle. 

II. La construction graphique de S n'est pas assez simple pour 
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permettre l'estimation rapide de l'erreur correspondant à une 
dislance h des ordonnées, quand on emploie la formule des tra- 
pèzes ou celle de Simpson. Mais il est facile de prouver que à est 
inférieur à Id et que, par suite, l'erreur, dans la première, est 
moindre que hd, et que | hd, dans la seconde. 

Pour le prouver, remarquons d'abord que l'on peut écrire ; 

tf = 2 [PQ — (QR — RS) — (ST — TF)], 

Or, on a PQ = d, et PQ > QR > RS > ST > TF. Pour établir 

ce dernier point, prouvons, par exemple, que PQ ou d est supé- 
rieur à d' ou QR. L'inégalité à établir, PQ > QR, devient succes- 
sivement : 

1111 






Soient, b, le poin! d'intersection de iB avee AC, A,, celui de ^K 
avec IL. On a 

-1 i 



Donc l'inégalité à vérifier devient : 

B6, + lik, > 0, 

et, sous cette forme, elle est évidente. 

Il est donc aussi facile d'estimer graphiquement une limite 
supérieure de l'erreur pour la formule des trapèzes ou pour celle 
de Simpson, que pour la formule de Poncelet ou pour celle de 
M. Parmentier. 
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Démonstration crai'hujue simultanée des résultats précédents. 
Pheuière extension au cas où le nombre des 

EST PAIR. 



H- Formule des trapèzes. On peut donner de la quantiié hd 
une interprétation graphique différente de celle du n" 7, appli- 
cable même iui cas où le nombre des ordonnées est pair. 




Considérons, par exemple, une aire2=«ABCDEF/' (Fig. 3) 
ayant encore sa concavité dirigée vers la base at, et ses ordon- 
nées aA, ôB,,..., F/, distantes entre elles d'une même quantité h. 
Cette aire est comprisejentre le polygone inscrit T = (iABCDEF/" 
et le polygone circonscrit aBjBgDiD^EiEa/"- 

La différence A entre ces deux polygones est plus petite que 
la somme des triangles 



BCBs, CDD,, DED,, 



EFEt. 



Par lo point B, menons BX, prolongement de CB, By, B5, Be, 
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Bés, Btp respectivement parallèles à DC, DiDD^, DK, EEs, EF. 
Les triangles énumérés tantôt sont égaux respectivement aux 
triangles 

ABB„ B,BX, yBJ, SBt, EjB)>, 

indiqués sur la figure par des hachures. Leur somme est infé- 
rieure au triangle ABtp, qui est égal lui-même à ia somme des 
deux triangles BPA, EQF, BP, EQ étant des parallèles à af. 
L'aire 2 est donc inférieure à l'aire T augmentée de ces deux 
triangles, autrement dit, on a 

«ABCDEF/< aires <oPBCDEQ/; 

résultat extrêmement simple et non encore signalé, que nous 
sachions (*). 

Si la eourbe était continuellement croissante (Fig, 4) ou 
décroissante, de A à F, la différence A est encore inférieure à 
un triangle ABif, dont le côté Bç est parallèle à la dernière 
corde du polygone T. Ce triangle ABcp est égal à la différence 
des triangles BPA, EQF. Par suite T -l- ABP— EQF, est une 
limite supérieure de 2, qui est encore égale au polygone 
aPBCDEQ/; comme dans le cas précédent. 

Pour une figure donnée, il est évidemmeni plus facile de 
construire les deux triangles ABP, EQF que le petit rectangle 
PQQ'P' de Poncelet (Fig. i). 

Dans le cas de la figure 5, on a d'ailleurs 

Dans le cas de la figure 4, 

APQ = 1 ft {y, - y,), DEF = i A {y„ - 1,„„.) . 



f) Si les ordonnées n'élaienl pas équïJialanles, il esl clair que l'erreuv conimise en 
prenant S=T est moindre qu'un triangle semblable à AB^ et ajai» pour hauteuf, 
au lieu de APt=Ji, la dislance das deuï ordonnées eonsÈcuLives de la courbe les plus 
iloignées l'une de l'autre. 
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Donc ABqj, somme des deux triangles dans le premier cas, 
différence dans le second, est égale à 

2 A (,'/» -!/>)-+- 2 '* (y«-< - !/J = /* [2 iy^ -^ y«-i) - ^ (j/i -+• .vJ ! - 

quantité que l'on peut désigner par /îDj 




Appelons T' Taire aPBCDEQf, de sorte que T' = T-i-/iD. 
On aura : 

T = /* (- ïi + ïi + J/î + ■ ■ ■ + y„-. "^ 2 ^7 ' 
T' = A (-i/î + Sï -H !/, -1- ■ • ■ + »/„_, -t- -y„_,)- 
On obtient donc une limite supérieure T' de l'aire 2, en 
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remplaçant dans l'aire du polygone inscrit T, la première el la 
dernière ordonnée par la seconde el l'avant-dernicre (*). 

1». Formules de Poncekt, de M. Parmentier et de Simpson. 
Ces formules n'existent que dans le cas où le nombre des ordon- 
nées est impair. Considérons donc seulement l'aire S=aABCDEe 
(Fig. 3). Le polygone m' inscrit de Poncelet est aABDEe. Pro- 
longeons DB jusqu'à sa rencontre en |3 avec «A. La différence 
M=M — m, entre le polygone circonscrit et le polygone inscrit 
est égale à la somme des quatre triangles 

ABB„ Bc.Bs = B,Bj3, D,DjC„ = jSB-f, DED^ = m, 

ou au triangle ABe. L'erreur maxima, dans la formule de Pon- 
celet, est donc moindre que \ ABs; dans celle de Parmentier, 
elle est inférieure à | ABî, et souvent à ^ ABe. 

On a vu, au n" 9, que dans la formule de Simpson, s < Aj3 cl 
/(5=2 (M — t). Dans le cas actuel, d'après le n° 11, M — t est la 
somme des triangles compris à l'intérieur de ABï et marqués 
par des hachures. Donc A(5<2 ABe, |fi(J<f ABe et, par consé- 
quent, l'erreur i est moindre que | kd. 

13. Première extension au cas d'un nombre pair d'ordonnées.. 
Ce qui précède suggère un moyen d'étendre les formules de 
Simpson, de Poncelet et de M. Parmentier au cas où le nombre 
des ordonnées est pair. Soit fi calculer approximativement l'aire 
S^S-f-CjO- désignant l'aire curviligne eEF/". On a 

( < s < f -1- ABe 
trapèze eEF/' = trapèze eEl-y< trapèze eEF/'-i- EE^F. 



(1 La formule des trapèzes (comme aussi celle de Simpson el les autres] est omplojSe 
à l'eslîmalion des volâmes, des ceotres de gravité, des momenls d'ioertie, etc. Dans le 
jaugeage des natires, par exemple, s,, y,,...., y» représentent les aires de aeclions âqui- 
distant«s et on peut les obtenir sur des coupes bien dessinées, au moyen do planimôlre. 
On conçoit aistaent qu'an moyen de deux iectures cel inatruraeni donne même immé- 
diatement T. D'aprÈa la tomiule démontrée dans le texte, on peut do même trouver T'. 
Donc, au mayen de qualre lectures, on peut obtenir une limite supérieure et une limite 
intérieure de la capacité du navire. 
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Ajoulons, en remarquant que 

trapèze eEF/"= - k {i/„^, -t- y„), ( + trapèze eEFf= T, 

AB« -I- EEsF ou ABe -t- îjBf < ABp ou Al"). 

!1 viendra : 

1 

T <s ->--fe (y.,_ , + )/„)< T-t- AD ou T'. 

Donc, en complétant la formule de Sinnpson, au moyen d'un 
trapèze inscrit, dans le cas d'un nombre pair d'ordonnées, l'er- 
reur est moindre que AD; autrement dit, la valeur trouvée, 
comme 2, est comprise entre les polygones aABCDEF/', 
ePBCDEQ/: 

Dans le cas des formules de Poneeiet et de M. Parmentier, 
on prouvera de même que T' surpasse p ou p' augmenté de 
~ h (^„_i+ J/Oî '"''is ces quantités peuvent être inférieures à T. 
Néanmoins, on établit encore aisément que les formules approxi- 
matives 

i i 

comportent une erreur -noindrc que AD. On a, en effet, 

1 
S = p di une partie de - ABe 

a^=- hly,,-! + y,) -i- une partie de c^B^. 

Donc, en ajoutant, 

1 
S^-[I = 2 = p-^- - fi (i/„_,-*-»;„)± une partie de ABf. 

De même, 

\ 
2 = p' H — h{y„^, H- y„] ± une partie de ABp. 
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Nouvelle limite supérieure de S. Formule nouvelle. Tboesiéme 

DÉMOnSïRATlON DE LA FORMULE DE SlMPSON. 

14. Nouvelle limite supérieure de S, I. Par les extréniilés 
G,E,G, I des ordonnées de rang impair autres que les extrêmes, 

menons, à la courbe (Fig. 5), les tangentes terminées en C, et C^, 




Ef et E^, G) et Gj, Ij et I^, aux ordonnées voisines. Prolongeons 
en outre CB, jusqu'à sa rencontre, en X, avec aA prolongé et IK 
jusqu'à sa rencontre, en Y, avec IL prolongé. On aura : 



aX4-r,G = 26B, lY ■*■ H =^ ikK, 



I, successivement, 



uX +^; = %, lY -i-y, = 2!/,„, 
«X = 2jr, — ys , ;Y = 2(/,a — i/,. 
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La somme M' des trapèzes 

0XB6. 6C,Csrf, rfE,Ej/; fGfiA kl,Uk, kKYl, 
est supérieure à S. Ces irapèzes ont pour mesure 

-h{Zy,-y,), 2%„ '2kij„ 2%,, .2%,, -ft(%,o-y,). 
Donc, après quelques réductions, 
M'=A(E+2I)+/i(»/,+t;,»-y,-y„)-/i|-(ys+yo--(y,-i-)/,o)|; 



ou encore. 



M' = A (E + 21) + Ihd — hd' = h^E-^- 21) -t- hd". 

II. On peut observer que d<Cà" <2rf. En effet, on 3 vu plus 
haut que d' est inférieur à d; donc rf"=2rf — d' est positif et 
supérieur à d. 

15. Formule nouvelle. I. On obtient une formule analogue à 
celle de Poncelet et où n'entrent que des ordonnées de rang 
impair, y^ et (/m exceptées, en posant 

S = p" [approximativement), 

p"=i(M' + m')=ME-^21) + ^Ad". 

On a, pour l'erreur maxima: 

,<_(M'-m') ou Ud-<_kd. 

On trouve, de même, une formule analogue à celle de M. Par- 
meniier : 

S =p"' (approximativemenl), 

p"'=i(2M' + m') = A(E + 21)-f-^Arf", 

2 4 

f < - kd" < - kd. 
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II. Nous verrons plus loin t^ue les deux foiiniiles précédentes 
ne soni pas très exactes. II n'en est pas de même de la suivante, 

S ;= iTi (approximativement), 

12 2 12 

^„ = p-" hd" = /* (E + 21) -1- - hd- kd", 

- ' 2« -f- 1 .3 ^ 5 2n -+- 1 5 

qui est nouvelle et semble avoir une valeur pratique et théorique 
analogue à la formule S=7ri, donnée plus haut. Malheureuse- 
ment la limite supérieure de Terreur est assez élevée, savoir : 

- hd" lid", 

5 2» + 1 3 

2»i étant le nombre des subdivisions de la base al. 

16. Troisième démonstration de la formule de Simpson. I.Des 
inégalités 

M > S, M' > S, 
on déduit 

2M + M' > 5S, - {m + M') > S. 

Or, 

i(^>M ^M')==-/({E-h2I-*- 4P) + -M"==s h--W. 



Done, 

M -t- - kd'- > S. 



II. De iiièmej des inégalités 

m < S, ( < S, 
oii lire 

Or, 

- (m + 2[) = - A (E ^- 21 + 41») ~ - lid = s 
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m. Puisque S est compris entre les quantités 



S = s (approximativement), 

avee une erreur moindre que i hd", qui est la plus grande des 
deux quantités ^ kd, i hd". 

IV. Nous avons donc la limite de l'erreur que comporte la 
formule de Simpson sous trois formes différentes : 1° i-C^hS; 
2°e< |M; 3" e<iM". 

La seconde est la limite la plus grande, mais la plus facile 
à construire, La première est au contraire la plus difiîeilc à con- 
struire, mais elle est toujours inférieure à la seconde, et, proba- 
blement aussi, à la troisième, car 

A<î=M — m', /trf"^=M' — m', 

et, en faisant un tracé graphique, on reconnaît que, le plus sou- 
vent M est inférieur à M'. Enfin, la troisième, certainement plus 
petite que la seconde, est presque aussi facile à construire (*}. 



Nouvelles démonstrations de la formule des tkapézes 
dans tous les cas. 



17. Cas d'un nombre impair d'ordonnées. I. Le procédé qui t 
servi au numéro précédent à trouver de nouvelles limites de S e 
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une seconde démonstralion de la formule de Simpson, peut évi- 
demment donner d"aulrcs limites et d'autres formules en grand 
nombre. Voici, par exemple, une nouvelle démonstration de la 
formule des trapèzes, c'est-à-dire une nouvelle estimation de 
l'erreur que l'on commet en l'employant. 
II. Des inégalités 

M > S, M' > S, 
on déduit 

M-)-M'>2S, -(M-i-M')>S. 
Or, 

- {M -1- M') = - fi (E + 21 + 2P) + - kd" = ( 4- - hd". 

On a donc : 

1 

( + - kd" > S > (. 

Par suite (4, 1), si l'on fait 

S = t [approximativement), 

on commet une erreur, par défaut, inférieure à i kd". 
IH. Posons encore : 



S = (' [approximativement). 

D'après le principe II, l'erreur maxima que comporte celte 
formule des trapèzes corrigée est inférieure h ^ hd" et le sens 
en est inconnu. 

18. Cas d'un nombre pair d'ordonnées. Prolongeons la 
courbe AL (Fig. 6), jusqu'à une douzième ordonnée AX=y,5, 
située à une distance h de Ll. Appelons a l'aire curviligne ILAi, 
et cherchons une valeur approximative de S-)-ff. Pour cela, 
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prolongcniis, la eordc KL jusqu'à sa rencontre, en ^, avec ).A. 
L'aire IL(). surpasse a. On aura, d'après le n" 14, 



(H» = -*(5j„-!,,.). 
Ajoulaiil eetle aire à M, qui est plus grande que S, il vient : 

h (ây, + %. + 2y, + ^y, + l S.0 + 1 y,) > s + -T. 

On trouve, de même, 

^'(^y> + l^^ + 2?/, + 2y, + 21/, -H 2y..) > S + <r. 

D'oii, en ajoutant et divisant par 2 : 

1 /7 ^ 3 7 \ 




Le premier membre de cette 
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Posons : 

1. = ft (- y. + Vî + î/î + ■ - ■ -t- !/,„ M- i/„ 4- - y„j, 

rf'; = 2rf, — d;. 
On pourra écrire la dernière inégalité : 
i 

ïi-4--m;> s -(-<t. 

La i^uantiié t^ est la mesure de la somme des trapèzes intérieurs 
inscrits dans S-\-<t. On a donc : 

(1 < S -«- ff. 

On déduit des inégalités précédentes, ces deux conséquences : 
1" Si l'on fait 

S -4- a = (1 (approximativement), 

on commet une erreur, par défatit, inférieure à \ hd[. 2" Si l'on 
fait 

S -t- CT = (, + - Mi, 



on commet une erreur de sens inconnu, inférieure à \ kd'[. 

l». Démonstration analytique de la formule des trapèzes 
dam tous les cas, sous la forme donnée au % V. I. Considérons 
Taire 2 d'une courbe dont la concavité est toujours tournée dans 
le même sens et dont on connaît n ordonnées éqnidistanies y,, 
y^t î/s.'-î !/»-»! y-, leur distance commune étant h. Posons: 

T = /» (2 2/i + ï/^ + ï' + ■ ■ ■ + .V~-^ -^ ^-' -t- 2 V')' 
\ \ \ \ 
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Les rcsullats des doux numéros précédenis penveni s'émrc : 

1 
T < 2 < T H- fiD — - AD' ; 

ce qui prouve que la formule des trapèzes 2=T (approximati- 
vement) comporte une erreur par défaiil moindre queftD^^ftD', 
II. On remarquera que D surpasse D', comme il est facile de 
le prouver (n" 1 0). On a donc, .à fortiori, 

T < 2 < T -t- AD, 

ce qui equndul au Ttsullat démon(ré graphiquement au § V. 
En \oiri une demonsliation analytique qui ne diffère pas essen- 
tiellement de celle du n" H (Fig. 1, à compléter par le lecteur). 
Ajoutons a un premiei trapèze analogue à aB,B^c, compris 
entre i/, et y^, des trapezts analogues à ^L^^ compris entre y^ 
c !/4î y* «•t i/jj, etc La somme de ces trapèzes sera : 

= 'n^p^ -^ y-^ + yi -^ ■ • ■ -t-y-^ + .y-» + g^-J =t + ad; 

ce qu'il fallait démontrer. 



Nouvelle extension des forbujles de Simpson, de Poncei.rt, etc., 

AU CAS d'un nombre PAIR d'ORDOHNÉES. FoRMCLE DE M. CatALAH. 

20. Estimation approchée de l'aire a. l. L'aire a étant plus 
grande que le trapèze ILaI (Fig. 6), on a 

1 

'^> â^(^<i + y<è' 
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D'après ie n" 18, a- est inférieur à IL^)., ou 

Toutes les valeurs intermédiaires entre les deux limites ici indi- 
quées peuvent servir à estimer approximativement a. Pour des 
raisons que nous indiquerons plus bas, la valeur la plus conve- 
nable est 

c = - (5/LAl -*- ILlX) (approximativement), 
G 

c'est-à-dire, toute réduction laite, 

L'erreur maxima est inférieure à 

^ {Ihly - (LA)) = 1 A (_ y,„ -H 2y„ - y.,), 

d'après le principe IV du n" 4. 

IL On peut trouver aisément une aire v inférieure au trapèze 
l\\Sk et, par suite, à j et telle que 

-(mil -t- ii) = -(5aAJ -1- ILU). 

On trouve sans peine : 

»-'.(|ï,.-».. + j!(,.)' 

On obtient géométriquement cette limite inférieure de a, limite 
au reste toute artificielle, en menant A^ parallèle à KL, Q^^ 
parallèle à Ka. Alors v=U^.^. 

81. Extension des formules de Simpson, de Poncelet, etc., au 
cns d'un nombre pair d'ordonnées. Formule de M. Catalan, 
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comme exleHsivn de celle de Simpson. I. Au lieu d'ajouler aux 
quantités s, p, p' {et aussi à ir,, p",p"', îia)» l'aire du trapèze 
ILa); pour avoir une valeur approchée de S-+-a, comme on l'a 
indiqué au n" 15, il vaut mieux y ajouter la valeur approximative, 

i 

— A(— ?/,o + 8)/„+ 5ija), 

trouvée ei-dessus. L'erreur commise dans l'estimation de S-f-cr 
est égale à la somme des erreurs commises en cherchant séparé- 
ment S et 0-. On prouve, sans peine, comme au n" 13, que cette 
erreur, dans le cas actuel, est inf'éneure, pour toutes les formules 
en question, à|/(D. 

II. La formule de Simpson devient ainsi ; 

1 

1 / 15 5 \ 

On trouve de même, pour une valeur approchée de la même 
aire : 

3 '' (4 y* -t- 52/2 + j «/s-*- %* ^- 4)/, + ■■ ■ . -f- 2y,o -+- 4t/„ + y„j- 

En prenant la moyenne de ces deux valeurs, on obtient l'expres- 
sion symétrique : 

^ ,/5 7 23 23 7 5 \ 

^A^- y, -1- i/s-i- Jî-t- yj-H !/„+ yc + f/, + f/3 + î/ii -^ y.-)-*- y.i -^- 2 ï'^) 
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La formule approximative, 

S = C, 

est due à M. Catalan. Nous allons la démontrer par une méthode 
au fond identique à la précédente, mais qui ne suppose pas qu« 
le nombre des ordonnées soit pair. 



Formule de M. Catalan, 

38, Formule de M. Catalan. I. Nous avons trouvé, aux n'" 1 '1 
et 19, que l'aire 

Ms = A (-)/,-(- y. + ■■■ +^„-. + y.-2-H^y-i)' 

dans le cas d'une courbe tournant en concavité vers l'axe al, 
supérieure à l'aire 2 de celte courbe, dont on donne les n 
ordonnées équidistanles j/j, y^,.—! V^- 

On peut regarder M^ comme la somme des trapèzes ayant 
pour aire 

dont le premier seul est un trapèze circonscrit de Poncelet, et 
dont les autres sont analogues à la limite supérieure de «r consi- 
dérée au n" 20, On peut aussi regarder M^, comme la somme 
des trapèzes ayant pour aire 
1 1 1 

II. D'après le n" 20, la somme des trapèzes qui ont poui- 
mesure 

/5 5 \ 1 

A(,vi-+-4/3), '^[i?h — ?/-'^iy>j 
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dont le premier est un trapèze iiiserii à 2, les autres sont ana- 
logues à la limite inférieure vde<T considérée au n" 20, est infé- 
rieure à 2. Celte somme est 

A (y, + ^ Ï2 + - ;/5 -H y* -H !/( + . . ■ + î/„_î -+■ - 1/„_, + - 1/„ j ■ 

Pour la même raison, les aires 

la 5 \ ■ ™ 

» T »-.-»-. -^i»-.' *(!/.->+!/.), I 






Inférieure à 2. 
Il en sera de même de la moyenne de ces deux sommes; 



«•.-'•\5Ï.+ 5!/.H-g!/.+S.-^-+!/.->-5Ï.-.-5!/.^.+ 5S.j' 

III. Nous avons dit, à propos de la formule de Simpson et au 
n° 20, et nous prouverons ultérieurement, qu'en général, la 
valeur la plus eonvenable à choisir pour chaque trapèze curvi- 
ligne composant 2 est égale auï: | du trapèze appartenant à la 
suite (1) ou (2), augmenté de ^ du trapèze correspondant de la 
série (5) ou (4). Donc, si l'on pose 

C=1(2M,-H»„), 

la valeur la plus convenable pour 2 sera donnée par la formule 
de M. Catalan, 

S^C (approximativement). 
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On irouve d'ailleurs ; 

M„ =T -t- /(D, m, = T — - hli tiD' , 

4 4 

= T-H~/(D — /iD'. 
4 

IV. Les principes généraux du n" 4 conduisent à une esii- 
niatioii de l'erreur beaucoup trop grande. On irouve £<C — m^, 

3 1 

C-~nh = -hU + -kO'; 
2 (1 

mais il est facile de trouver une limite moins élevée. Puisque 
D>D', on a C et 2, compris d'après le n° 19, enlre les quantités 



Donc l'erreur commise, en posant 2^=C, est moindre que la plus 
grande des quanliics 



-hù', T + /iD — -/iD'- 



e'est-à-dire évidemment la seconde. En pratique, on peut dire 
que s est moindre que ^ liD, 



SnCONDIJ POlillULE DE SlMPSON. 

S3. Élément de la seconde formule de Simpson. I. Considé- 
rons la parlie abcdDCBA. de l'aire S, et enfermons-la entre deux 
limites. 

Cette aire est plus petite que la somme des trapèzes 

«tBX, bC,C,d, 
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dont il est question an ii" 14. Laire de ceux-ci est 



On peut remarquer que le second de ces trapèzes a h même 
aire que le irapèze 6BD'<?, D' étiint le point d'intersection de BC 
avec Drf. On a donc : 



résultat qu'il est facile d'élàblir directement. 




II. L'aire courbe «fierfDCBA est supérieure à l'aire polygo- 
nale aècrfDCBA, qui a pour mesure 

i 

- A (y, + 2!/i -1- 2j/, -V- !/0- 

III. Nous verrons ultérieurement que la valeur la plus con- 
venable pour l'aire courbe oABCDrf est 

in 5 "I 3 
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L'erreur commise esl plus petite que la plus jurande des deux 
différenees, 

- 4(3}, + 3y.) - - ;. [ ji * 3y, + j. + j,] 

, r 1 1 1 in 

j /> (j, H- 3j, + 5j, -4- 5.) - - ft (j, + %, + 2j, H- J,) 
-'•[-5Ï.-5S.+ i!/.-j!/.]' 

c'est-à-dire r^ue la première. Celle quamîlé esl facile à inter- 
préter géoniéiriqnemont. 

34. Seconde formule de Simpson. I. Supposons de même que 
l'on ait a chercher l'aire 2 d'une courbe dont la concavité est 
tournée vers l'axe des a; et dont on connait (5n -+- 1 ) ordonnées 
équidistanles yj, y^, «/j,..., y,„+,. On a, comme limite supérieure, 

-A(2ys-H 2i/i, -t- 2i/6-*- 2j/c h- 2^8 ■•- ^o + ■■■ + 2)/,„_, -+- 2i/,„), 

et, comme limite inférieure, 

1 

^ /( (y, + 2i/2 + 2»/3 -1- ■ ■ ■ + 2y5„ H- i/^^,). 

II, On prendra pour valeur approximative de l'aire, le quart 

de la première expression plus les trois quarts de la seconde, 
c'est-à-dire ; 

La formule 

8=:^' (approximativement) 

esl la seconde formelle de Simpson. L'erreur maxima commise 



y Google 



est plus petite que la différence de s' avec la limite supérieure. 
On a donc 

K|ft[-,Vi+y3+V=-%<-*-?/»-^ïo-2y'-+---%=-î-t-î/3.-.+?/a.-.V-+i]- 

III. On peut obtenir une autre limite de l'erreur, plus facile 
à calculer à priori, par le procédé qui a servi plus haut pour 
la formule des trapèzes, au § V. On trouve immédiatement 
e<|ftD, D étant donné par la formule 

1 i 

OU représenté par ia somme ou la différence de deux triangles 
analogues à ABP, EQFde la figure 3. 



Formule de Weddle. 

S5. Limite supérieure et limite inférieure. Considérons une 
courbe ABCDEFG concave vers l'axe ag, et sept ordonnées Aa, 

B6, G^r, équidistantes d'une quantité h, que nous désignons 

encore par 1/1,1/3,..., y^. Estimons l'aire Sj^^aABCDEFGg'. 

Par les extrémités B, D, F des ordonnées de rang pair, me- 
nons encore les tangentes B1B2, DiDj, FiF^ limitées aux ordon- 
nées voisines. La somme M| des trapèzes 

«B,BsC, cDJ)ie, eF,F,3, 

est supérieure à S,. Or, 

Menons les cordes AC, CE, EG, Le polygone «ACEG^ est 
inférieur à Sj et a pour mesure 

m, = A(j/, -H 2»/5 + 2i/, + !/,}■ 
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Menons aussi les cordes AD, DG. Le trapèze «ADGy esl aussi 

inférieur à S] et a pour mesure 

w*3=-/i(!/.-v-2.y.-i-î(,). 

Des inégalités 

m,<S,, mi<S,. 
on déduit 

i 

jim, -H âflis < 3S, , - (3»(, -+- âm^) < S,. 

Or, 

1 5 

,H,= - (3m. -.- im,) = ^ /i (2(/, + 2)/5 + 2»/* + 2y, + 2f/,). 

8e. Formule de Weddle. Laire Si élant comprise entre Mj 
et ïWj, on peut prendre pour valeur approximative de celle aire 




toute valeur comprise entre M^ et mj. La valeur la plus conve- 
nable est, comme on le verra plus loin, 

tv= - {3M.-1- m,)=~h{y, + %, + y, + Ci/, -h y. + %, + y.). 
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Nous poserons donc 

S, = w (approximativement), 

ce qui est la formule de Weddk. 

L'erreur commise est plus petite que 



-Am,^m.)- 



- 3 (M, - 



On peut déduire de ce qui précède une formule servant à 
trouver les aires dont la base est divisée en &n parties égales. 



Exemples ncsiériques (*). 

ST. Premier exemple. Soit à chercher Taire de la courbe dont 
l'ordonnée esl égale au logarithme népérien de x, diminué de 3, 
X variant de lOOOà lOOttetA étant égal à l'unité. On trouve, en 
dix-billionièmes (l'aire exacte étant donnée par le calcul intégral) 



»• 


— 13 407 74,79, 


V 


= se 77â1b.51, 


ïl 


:= 430 093 30.20, 


V 


= 175 57128.09, 


•/ 


= 216006 17.S7, 


» 


= 259 79807.20. 




d =10793, 




d'= 647», 




rf"= 15110, 




1 =12951. 



= 780)7128. 




= S, .... 


i < 4317 


= S, .... 


e < 647S 


= s, . . . . 


E < 129S1 


= S, .... 


e < 4396 


= S — 2158 . . 


s < 6476 


-S- 1080, . 


£ < 5396 


= S -1- 720, . . 


e < 7195 


==S 


E < 6476 


= S-1079, . 


, < 7SS5 


= S -t U39, . 


E < 10073 


= S 


£ < 8654 



(') Les exemples traités ici ont exigé de nombreux calculs, pour le premier i douze 
décimales, poar les autres à huit déciuiales cl avec deux ou trois groupes de données dis- 
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Les calculs ont été faits en allant, jusqu'aux inllionièmcs, c'est-à- 
dire, en prenant deux décimales de plus que celles qui sont con- 
servées. La première formule de Simpson et celle de Weddie 
donnent un résultat exact jusqu'aux trillionièmes. Même si l'on 
n'ennploie que trois ordonnées, pour la première formule de 
Simpson, c'est-à-dire si l'on pose s=!/, ■+■ 4j/4 ■+- y^, la différence 
S — s est moindre qu'un cent-billionième. 

Quand les erreurs théoriques peuvent être calculées de plu- 
sieurs manières, nous n'avons inscrit dans la dernière colonne 
que l'erreur maxima la plus faible. 

28. Autres exemples. Dans le cas précédent, p, p', tt^ sont 
calculés au moyen de cinq ordonnées seulement ()/5 et y^ étant 
exclus), p", p'", itj au moyen de 6 (j/4 étant exclu); les autres 
valeurs ont été trouvées au moyen des sept ordonnées. Voici 
trois autres exemples, où, pour chaque formule, on a employé 
seulement sept ordonnées ; h = \ pour les S premières, h^-^ 
pour les trois suivantes ; enfin A=| pour les dernières. Tous les 
résultats sont donnés en millionièmes. 







Coiirbe 1. 


Couitc II. 




Courbe ]!l. 


Aire exacte . . 


S =6i)5147 


S =78S598 


S 


=19152^3 


Simpson 1 . 




S =S+25 


s =S 


S 


= S-5 


Simpson II. 




s' =S-t-48 


s' =S-2 


s' 


==S-1 


Weddle . 




w =S + ^ 


IV =S-1 


ÎU 


= S 


Catalan. . 




C =S-.-fi6 


C =S-43 


c 


= S-10 


Trapèzes . 




( =S + 1730 


( =S-H57 


( 


= S — 134 


Poncelet . 




p =S-+-376 


p =5-251 


P 


= 8—327 


Parraenlier 




p' =8-163 


p' =S-t-124 


P' 


= S-Hl02 


Dupain, . 




TT, c=S-43 


TT. =S+45 




=S-t-7 




jî"=S + 603 


p"=S-ZQ\ 


P' 


= 8-468 




;,"'==.S — 488 


p'"=S-+-387 


p' 


'=S-t-377 


Formule nouvelle 


^, =S — 5 


jii =S+S5 




= S-v-2 


tinctcs. Il se pourrait dan 


qu'il se fût glissé quelques erreurs dan 


'Z 


vésukats.inais 


elles ne sont i)as bien cou 


dérablas,pensons-nou 


. La eourbe 1 s sa convexité (ournéo vers 


Va\e al, aalremea 


diî,rf, 


d', ci", J sont négatifs 
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La courbe I est une hyperbole équîlatère, de sorte que l'or- 
donnée est égale à [1 : ce], et x varie de 1 à 2. Dans la courbe II, 
rordoiinée est égale à [1 : (1 -i- a;^)], x variant depuis jusque 1 . 
La courbe III est une circonférence de rayon égal à 2; on es- 
time l'aire comprise entre cette circonférence, un premier rayon, 
un second rayon perpendiculaire au premier et une parallèle au 
premier par le milieu du second. On peut prouver, par la géo- 
métrie élémentaire, que les aires de ces trois courbes sont pour I, 
le logarithme népérien de 2, pour II, f ît, pour III, s^^ + î 1/3. 

S». Exemples dans le cas d'un nombre pair d'ordonnées. 
Voici, pour une partie des formules précédentes, les résultats 
que l'on trouve pour les courbes I el II, en se servant de huit 
ordonnées. Pour les formules de Simpson et de Weddie, on les 
applique à l'aire comprise entre la première et la septième or- 
donnée; on ajoute ensuite au résultat irouvé une valeur de u 
calculée comme il est dit au n" 20. 





tourbe 1. 


Courbe 11, 


Aii-e exacte . 


. S =693147 


S =789398 


Simpson 1 , . 


. s = S-+-4 


S =S-h4 


Simpson 11. . 


. .5'=S-t-18 


s'=S-v- 16 


Weddic. . . 


. w=S-8 


«;.= S+ 13 


Catalan . . . 


. C = S + 59 


C = S — 22 


Trapèzes . . 


. T = S-t- 1-272 


( = S — 8S0 



Choix d'une valeur ja intermédiaire entre une limite supérieure 
et une limite infëbieube de s, comme valeur approchée de s. 

30. Principe. Dans les paragraphes précédents, nous avons 
enfermé Taxe S à mesurer, entre celles de divers polygones 
M, M', elc, plus grands que S, et d'autres polygones plus petits 
m, m', elc. Appelons, en général, L la limile supérieure de S, l la 
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limite iid'érieure. Ces polygones L, l sont des sommes de tra- 
pèzes ayant pour hauteur commune une longueur k, aussi petite 
qu'on le veut. Les différences L — S, S — l, positives pour les 
eourbes à concavité tournée vers la base de l'aire à mesurer, 
négatives dans le cas contraire, ont l'une et l'autre ponr limite 0, 
quand h décroît indéfiniment. 

Soit R, le rapport inconnu, toujours positif des différences 
L — S, S — /. De la relation 



-R, 



Lorsque h décroît, en général, le rapport R varie. Dans beau- 
coup de cas, on peut en déterminer la limite r, pour h=0, sans 
connaître la valeur exacte de S. Supposons que nous prenions, 
pour valeur approchée de S, h quantité p, définie par la relation 



L'erreur réelle e = f^ — S sera donnée par la formule 
_L-^-rl L + Rl _{l — l){l\ — r] 

Le dénoniinateur de ■; est supérieur à l'unité; le numérateur 
est le produit de deux facteurs dont chacun tend vers zéro en 
même temps que h, puisque lim L =lim/=S, limR = r. 
Si, dans l'expression de (x, on mettait au lieu de r une autre quan- 
tité positive, l'erreur aurait pour numérateur le produit de 
deux quantités dont une seule, L — I, aurait pour limite zéro. 

Donc, dans les cas dont il s'agit, pour h suffisamment petit, 
il est préférable de prendre pour S la valeur jx donnée plus haut. 
Dans les autres cas d'ailleurs, en général, on n'a aucune raison 
■pour ne pas faire de même. 
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L'erreur e, ramiiie on l'a dil au n" 4, IV, csi plus petite que 
la plus graixie des deux différences 



Si r > I, s < L — p; sir < 1, î < fj, — (. 

31. Cas où le principe précédent est applicable. II existe un 
cas assez général, eomme on le verra au § suivanl, où l'on peut 
trouver r aisément et voir que p. est la valeur approchée de S la 
plus convenable. C'est celui où la courbe considérée est continue, 
où l'inclinaison de ses tangentes et sa courbure varient d'une 
manière continue et où, enfin, ta courbure ne varie guère, soit 
pour la courbe entière, soit pour les diverses, parties dans les- 
quelles on est amené à la décomposer. 

Analjtiquement, si y.= F'x est l'équation de la courbe, 
F'x, F"x, F'"x doivent être des fonctions continues et F"'a; varier 
peu entre les valeurs extrêmes de x que l'on considère. 

On observera que souvent l'on pourra vérifier rapidement ces 
conditions, soit au moyen d'un tracé graphique, soit au moyen 
de calculs assez simples, quand F'x ne sera pas une fonction 
trop compliquée. 

Dans le cas où l'on ne parvient pas à voir si ces conditions 
sont remplies, et où aucune propriété spéciale de la courbe 
ne force à choisir pour (x une valeur autre que celle qui est 
donnée plus haut, le mieux évidemment est de s'en tenir à cette 
valeur, qui a, sur toute autre, l'avantage de se justifier dans le cas 
général indiqué ci-dessus. 
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CALCLIL DE 7'. 

33. Cas de trapèzes simples. 1. Considérons d'abord le tra- 
pèze curviligne aAB6 (Fig. 9), Appelons la première ordonncie 
F'{a;), la dernière F'(x + h). Ce trapèze a pour mesure 

p(x + A) — F{x) = /-F' + J- A'F"^ i AT "(a;,), 

F', F' éiant écrits, pour abréger, à la place de F'(x), F"(a;), et a-) 
étant une valeur intermédiaire entre x el a; -4- h. 




Par le point A menons la tangente AA^ terminée en Aj à B6 
prolongé. On aura 

ùAs = aA + AF"(j:] = F' + /(F". 

Par suite, le trapèze aXA^b, circonscrit au trapèze curviligne 
«AB6, a pour mesun; 

1 A[F' 4- (F' -+- /iF")] = h¥' + ^ h'F-. 
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Le trapèze inscrit aABb a pour mesure 

- h [F'(x) -h F'[x -H /()] -= kr ■+■ - k^F" -h - h?V"' (x^), 



x^ étant compris entre œ et a: - 
Dans le cas actuei, 



-F'UxA 
e. ^ " 



Par suite, si A tend vers zéro, et, par suite a:) el x^ vers se, on 
a lim R= r= 2. De plus, si F"'(a;) varie peu de a; à a; + h, 
R a aussi une valeur voisine de 2. 

On trouve le même résultat si Ton remplace le trapèze cir- 
conscrit «AAa6 par le trapèze aB,Bft, BjB étant la tangente 
en B, prolongée jusqu'à sa rencontre avec l'ordonnée «A. 

II. Soit (<i le point d'intersection de la corde AC avec 6B. 
On aura 

6&, -i(«AH-cC} = i[rM-i-F(^:^2A)] 

= i[2F' + HhV" + 2A'F'"(x5)]= F' + hV -t- fe'F'"(i,), 

Xs étant compris entre x cl x -\- 2A. L'aire du trapèze akh^b 
sera 

Si l'on prend celte aire poin' limite inférieure, dans la valeur 
de R, il vient : 

G ^ " 2 ^ ' 
valeur peu différente de ^, si F"'(x) varie peu de x à x -i- 'ih . 
Si k tend vers 0, on trouve dm R ^ ?■ = i. 

On arrive au même résultat encore, si l'on remplace le tra- 
pèze circonscrit aAA^ô par aBiBb. 
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m. Soil X h point d'inlersectioii de CB avec Aa. On a, 
comme on l'a vu an n° 14, 

qX = 26B — cG = W(x -t- h) — F'(x -4- 2A). 
Le trapèze nXBba pour mesure 

iA[5F'(x-t-A)~F'(x4-2ft)] 

= 1 A [3F' + 3ftF"-t- - fe'F"'(a;,) — F' — 2AF" — 2fe'F"' (x»)] 

= AF' -t- 1 A'F" + j A'F"'{x,)— A=F'" (x^), 

X4 étant compris entre a; et a; 4- A, % entre x et a: H- 2A. 

Si l'on prend celte aire pour limite supérieure, dans la valeur 
de R, et a\Bb pour limite inférieure, il vient 



L — S 



jF"'(x,)-F"'(x»)-iF"'(xO 



1f-"(xO-^f"'(.,) 

valeur voisine de S, si ¥'"(3:) varie peu de a; à a: + 2A. Si A tend 
vers 0, on trouve Hm R = r ^=5. 

33. Somme de trapèzes. I, Dans l'élément de la seconde for- 
mule (le Simpson, faisons y, = F' (x), y^ = F' (x -l- A), etc. 
L'aire de cet élément sera 

F(x + 5A) — Fa; = 3/tF' + - A'F"-+- -A"F"' [x,), 

Xf, étant compris entre x et x 4- SA. 
La limite supérieure est 

-A[F-(;c-i- A) + F'(x+2ft)] 
= I &[2F' + 3fcF" -t- A'F'"(x,) + 2A'F"' (a;,)] 
= 5AF'-i- - A'F" + - A'F"'(x,) + 5ft'F"'(x,), 
x-j étant compris entre x et x 4- A, Xg entre x et x + 2A. 
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La limite inférieure esi 

- h[F'3: + 2F'[a: + h] + 2F' (ï -t- 2A) -t- F'(x -t- 3A)] 

= 3feF'-t-^ft'F"-+-^ft'F'"(a:,)-v 2ft=F'" (x.„) + ^ /i'F'"(x„), 

iCg étant compris entre a; et a; + A, xk, entre a; et « + 2A, 3:,i 
entre X et a: + 3A. 
Par suite, 

R=^ ^ , 

- F-"{x,] - - F"{x,) - W M - - F"'(a:„) 

valeur qui a pour limite r =^ 5, et qui est d'ailleurs voisine de 
5, si F'" X varie peu liex t x -i- 3ft, 

II. Des calculs absolument analogues prouvent que, dans le 
cas de l'élément Weddle, r = ^ ei que R a une valeur voisine 
de j, pourvu que F'" vm-ie peu de \ àx -t~ 6ii. 



Application aux diverses formules précédentes. 

34. Formules de Simpson, de Weddte, et de M. Catalan. La 
première valeur de r trouvée au n" 32 prouve que la formule de 
Simpson, sous la forme 

S = s (approximativement), 



M ayant le sens indiqué au n° 6, est bien établie d'après le prin- 
cipe du § XIII. On voit de plus qu'elle sera probablement très 
exacte, si F'" (x) varie peu dans chaque intervalle k. 
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La seconde valeur de r (n" 32) conduit i 
sous la forme 



- = -(2M + w*'), 



m' ayant le même sens qu'au n° 8. On ne voit pas aussi bien, au 
moyen de celte valeur de r, pourquoi la formule de Simpson 
est généralemeni si exacte; car, pour établir que R est voisin de 
r = ^, on doit supposer ici que F'" (x) varie peu dans chaque 
intervalle 2/t. 

La formule des trapèzes n'est pas aussi rationnelle que la for- 
mule de Simpson, puisqu'elle suppose que R, dans ce second cas, 
est égal à 1, valeur qui n'est pas voisine de |. 

La troisième valeur de r, combinée avec les considérations du 
n" 20, II, conduit naturellement à la formule de M. Catalan; la 
quatrième à la seconde formule de Simpson, la cinquième h 
celle de Weddle. On admet que F'" varie peu, dans un inter- 
valle 2A pour la première, SA pour la deuxième, 6ft pour la 
troisième. 

35. Formule de Poncelet, de Parmentier, de Dupain, etc. 
l. Pour une aire divisée en 2w trapèzes curvilignes, les limites 
M et m données au n" 6 sont telles que dans le premier et le 
dernier de ces trapèzes, la valeur de R est voisine de 2, pourvu 
que F'" (x) varie peu dans un intervalle h; la valeur de R est 
voisine de j, pour les (2h — 2) autres trapèzes curvilignes, 
pourvu que F'" (x) varie peu dans un intervalle 2ft. 

On reconnaît aisément qu'il n'y a pas moyen d'employer la 
valeur 2 pour estimer le premier et le dernier irapète, et la 
valeur i pour les autres, parce que l'on introduirait ainsi dans 
la formule les dérivées des ordonnées extrêmes. Poneelet a choisi 
une valeur intermédiaire, savoir i, en posant 



1 

p = - (M -I- m). 
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M. Parmenlier, remarquant que le nombre des trapèzes inté- 
rieurs esl en général supérieur à celui des extérieurs, a fait r = ^, 
ce qui lui a donné 



l 



-=-(2M + m). 



\ - 



II. Nous sommes arrivés.à la valeur tc^ , plus rationnelle encore, 
de la manière suivante. Appelons u, la petite aire analogue à 
ANB comprise entre un arc ANB de la courbe, et sa corde AB 
dans un intervalle h (Fig. 9); l'aire ANBA3, comprise entre la 
tangente à une extrémité de l'arc et la courbe, sera sensiblement 
2«, puisque la limite de leur rapport est égal à 2. L'aire analogue 
à ANBÔ] entre la courbe ANB et la demi-eorde Aôj correspon- 
dant à un intervalle 2A, sera, pour la même raison, à peu près 
égale à 4m. On aura donc approximativement 

M = S-t-2n.2«=.S -t- 4rtM, 
m = S — M — {2w — 2) 4m — M = S — [8m — 6)« , 

On tire de là, en éliminant tt, les relations approximatives 

(4n — 3)M -t- 2»»t. = (6w — 5) S, 
1 \ 



{in — 3) M M- 9«m 



2m- 



1 posant, comme au n' 7, 



TT, = - (2M + m) (M — î«)=M hd hd. 

Les valeurs 2m, Au, données plus haut, représentent exacte- 
ment les aires analogues à ANB, ANB?*,, dans le cas d'une para- 
bole du second degré. Il en résulte que la formule S =ît, donne 
exactement l'aire d'une parabole du second degré, et, par suite, 
d'après un principe général que nous établirons dans le para- 
graphe XVII, aussi celle d'une parabole du troisième degré. 
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m. Des considérations analogues à ct'iles qui précèdeni prou- 
vent que les valeurs p", p'" du n" iS peuvent èlm remplacées 
par la valeur plus rationnelle . 

= /( (E + 21) -H - Arf" '■ hd". 

La formule nouvelle, 8 = %, est exacte aussi dans le cas d'ui 
parabole du deuxième ou du troisième degré. 



Calcul df. l'rrredr par le théorème de Tatlor. Foriujles 

DE SlMPSOX ET DE WeDDLE. 



36. Hypothèses préliminaires. Dans ce qui suit, nous allons 
ajouter les hypothèses suivantes à celles que nous avons faites 
jusqu'à présent; J° Les dérivées F"x, F'cr, etc., qui vont entrer 
dans nos calculs ultérieurs, sont continues pour les valeurs de x 
considérées. 2° La dérivée d'ordre le plus élevé parmi celles que 
nous emploierons varie peu pour les valeurs de ce considérées. 

L'erreur commise en employant les diverses formules données 
plus haut est de la forme 

E = aP-{a;] -*- b^i^:,) + cT'ix,) -+■ etc., 

F'x étant la dérivée d'ordre le plus élevé dont nous ayons besoin, 
ce,, aîg, iKj, ... des valeurs de x inconnues, mais comprises, entre 
les valeurs extrêmes de x. Cela posé, soit Xy, une valeur de x, 
comprise aussi entre ces valeurs extrêmes deic, on aura 
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Par hypothèse, la partie de e écrite sur la seconde ligne est 
très petite. On pourra donc écrire approximativement 

, = („ + 6 + c*...)F'(x,). 

Cela revient à mettre, dès le début, dans l'expression de s, ou 
même dans les quantités qui ont servi à trouver s, Xy, au lieu de 
x.],x^, ics, etc. C'est ce que nous ferons désormais, de sorte que 
toute égalité où entre x n'est qu'approximative. 

37. Formules de Simpson et de Weddle. I. L'élément de la 
première formule de Simpson, e'est-à-dire {h (y^ + %2 4- j/^), 
devient, en posant y^^F (ce), 

;-/*[4F'(x)-f. F'(x + /()-t-F'(a:— &}], 

et, en développant, 

3 3b 

L'aire exacte est 

F(a;-4-/i) — F(x — A) = 2feP"-t- -ft=F'"+ — F'(x^). 

L'erreur est donc approximativement 



= — h'F-(xu_) /t'F'{x„)= — A'F'fXu). 



Si F' (x) est du second ou du troisième degré, P (x) est nulle. 
Donc la formule de Simpson est exacte pour une parabole du 
deuxième ou du troisième degré. 

Si la base de l'aire à mesurer est 2/ift, l'erreur pour chaque 
élément de base 2A est de la forme ^ k^ F° (a:^ ) et ces erreurs 
s'ajoutent, puisque F' (x) variant peu, ne peut pas, en général, 
changer de signe. 
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II. L'élément de la seconde formule de Simpson peut 
s'écrire 

^&[F'(x) -- 3F> + A] -H W{x -,- 2k) + F'{x ^ 5fi)l 
Q Q 27 '^T 

2 2 8 16 ^ ^ 

L'aire esacie est 

9 9 27 81 

F(3:-t-3/() — Fx=5AF'+ - AT"+ ~ ft'F"'+ - A'F"+ — A'Pf»;.). 
^ ' 2 2 8 40 ^ '^ 

L'erreur est donc, approximativement, 

,= 5! ;(»F'{x„)— — A'F'(x„) = — A''F'('a:„). 
16 ^ ' ' 40 ^ '^ 80 ' ■^' 

Celte erreur est nulle dans le cas d'ime parabole du deuxième 
ou du troisième degré, 

III. L'éiémenl de la formule de Weddie est : 

— A [F'(x — ôh) +■ bF' (a; — '2h) -t- F' (« — A) -i- 6F' [x] 

-V- ¥'{x 4- A) + W{x + 2A) -)- F'(x + SA)] 
= 6F'+-9A'F"'-4- — A''F'-i--A'F'"(3-^). 

L'aire exacte est 



valeur nulle pour les paraboles du second,dulrorsième, du qua- 
trième et du cinquième degré. 
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IV. La valeur donnée par la seconde formule de Simpson 
a été établie de manière qu'elle fût exacte pour une parabole 
du troisième degré par Newton lui-même (Gauss, Werke, III, 
p. 202). 

La méthode qui a servi à trouver la formule de Weddie revient 
aussi à faire en sorte qu'elle donne l'aire exacte d'une parabole 
du S' degré dont on connaît sept ordonnées (voir plus bas, les 
notions historiques sur les diverses formules). Le second et le 
troisième résultat donnés plus haut, n'ont donc rien que de très 
naturel. 

Il n'en est pas de même du premier. La valeur approchée 
ï ft(yi -i- %2 + î/a) de l'éiément a été obtenue, comme nous le 
dirons dans l'historique, en cherchant l'expression qui donne 
l'aire exacte d'une parabole du second degré. Cependant, cette 
formule est exacte aussi, comme nous venons de le voir, pour 
une parabole du Iroisièmedegré. C'est là un théorème curieux (*) 
dont la formule de Gauss exprime, au fond, une belle générali- 
sation pour les paraboles d'ordre supérieur à ordonnées non 
équidislantes (Gauss, ÏFtr/tp, III, pp. 163-196,202-206). 

Il en existe une autre généralisation pour les paraboles quel- 
conques d'ordre pair, comme nous allons le montrer, avant de 
continuer la recherche de l'erreur pour les autres formules. 



XTII 

Théorèmes sur les paraboles. 
38. Théorésie I. Soit 

h (ay, -^ 6)/i -t- cj/, H + gy„), 

une formule donnant l'aire exacte d'ime parabole d'ordre p, p 



l'I Ce tbéorème, donné par M. Calalaii, dans les nouvelles Aniiatea de Mathémaiiquea, 
i" série, l. XVI, p. 312, et dans le Itanuel des Candidals à l'école polylechnique, l. Il, p. 
99Ë, a été recrouvépar M.P&rmentier (Congrès de Monipellier de l'Association française 
m des Sciences, iSlS). 
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étant inférieur à n, au moyen de a ordonnées de cette parabole , 
distantes entre elles d'une quantité h. Si on applique cette formule 
à une courbe y = F's, l'erreur commise s sera de la forme 

En effet, considérons une parabole d'ordre p, y = f (^), 
osculatrice à j/ 1= F' («), en un point quelconque [X, F'(X)], 
de sorte que 

/■' (X) = F' (X), / " (X) = F" (X), . . . , r+' (X) ^ V^i (X). (1 ) 

Soient 

a:, = X + kh, a;2 = a:, -t- /(, x-.-^x^ -t- k, ... «„ = a:„_, + /(, 

ies abcisses correspondant à )/„ ?/a, y^, —, y„- Posons : 

S = F (^,) - F (x,), 

S' = /i[aF'(x,)-H6F'(x.)-i---- -t- 9F'(a:,.)], 
s'=/i[af(x,)-,-b/-'(x,)+--- +sfrM. 

Dcveloppons, suivant les puissances de h, la différence S ■ — s, 
en ne laissant dans le développement d'autre abcisse que X. A 
cause des relations (1), on aura évidemment 

S~s = A,/('^'F'^'(a;p), 

car les autres dérivées de F, /'sont égales et f'^'' (x) = 0, Pour 
la même raison, 

S' — s' = AsA''+'F'^(a;^), 

Mais, d'après l'hypothèse faite sur la formule qui donne s', on a 
s'=s. Donc : 

E = S — S' = (A, — A4) h'-*^P'+\x^) = AA''+'F'^^(it^}; 

ce qu'il fallait démontrer. 
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En particulier, si p = («— 1), on a S — S' = A, A""*"' F (x^), 
ce qui est la justification de la méthode de quadrature approchée, 
due à Newton, dont ii est parlé au n° 1, pour les courbes dont 
l'ordonnée F' (x) vérifie les conditions énoncées au n' 36 (*). 

39. Théorème II. Soit 

A[a(y, -+-!/„) -t- 6 (^2 + y„-.) + e(y, + y„..,) -f- ■■■ ], 

Mlle expression donnant l'aire exacte d'une parabole d'ordre 2q , 
2q étant inférieur à n, au moyen de n ordonnées de celte para- 
. bole, cette expression étant symétrique par rapport aux ordonnées 
égalemetil éloignées des extrêmes. Si on applique cette formule à 
une parabole d'ordre 2q -1- 1, elle en donnera aussi l'aire exacte. 
Soient, en effet, y = F' (x), l'équation de cette parabole, et 
2X ^ a;] H- a;,,. On aura 

valeur qui, développée, ne contient que des puissances impaires 
de h. Ensuite : 

(.rF'(x-i(»-4)).) + F(x*l(„-l)4)l+ 1 

(' ["' ("^i*"-^''') * "' (^ * 5(«- 3)*)] *elc.| 

valeur qui, développée, ne contient aussi que des puissances im- 
paires de h. D'après le n" précédent, le premier terme de la dif- 
férence S = S' est 

ou plutôt. 



{') Si l'on veut justifier la mëlhode àe Sewlaa, en einplojant des tracés graphiques, on 
reconnaît Lien vile qu'à pariir du troisième degrÉ, très souvent les paraboles a'écaiient 
beaucoup des courbes avec lesquelles elles ont autant de luints communs qae possible. 
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[iiiisquc F' (x) est d'ordre (^q ■+■ 1) et F^"^ x constant. D'autre 
part, celte différence ne contient que des puissances impaires 
de k; donc A = 0, c'est-à-dire que S = S'. 
CoROLLAiHE. Si l'oii applique la formule 

à une courbe y = F' (x), l'erreur sera de la forme 

XVIH 

Calcul De l'erreur. Formule des trapèzes, de Poncelet, etc. 

40, Formules des trapèzes, de Poncelet, de M. Parmentier , 
de M. Dupain et formule nouvelle. I. D'après ce que l'on a vu au 
n" 32, l'erreur que comporte la formule des trapèzes est égale, 
approximativement, à 

fi — i] W'(x.)= ~/('F"' (x.). 

II. Dans la formule de Poncelet, l'erreur, pour les élémcnls 
extrêmes, est 

Ifi s! ^ ■' 24 ^ ■ ' 

Pour les éléments intérieurs, on a, d'après le même numéro, 



111. Dans la formule de M, Parmentier, l'erreur, pour les élé- 
ments extrêmes, est 
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Pour les éléments extérieurs, l'erreur, au contraire, s'exprime 
en A*, confime il est évident, d'après 32, 11 et 36, Si l'on prenait 
Terreur totale pour tous les éléments, on pourrait faire en sorte 
que l'erreur ne contienne que des puissances impaires de h, 
comme on l'a vu au n° 39. L'erreur dans la formule de M. Par- 
mcntior est donc du troisième degré en h, mais cette erreur en h'' 
provient uniquement des éléments extrêmes. , 

IV. Enfin les formules S = :t,, S = % étant exactes pour des 
paraboles du deuxième, et, par suite (n''39),du troisième degré, 
l'expression de l'erreur est de la forme 

41. Formule de M. Catalan. L'élément t dont il est parlé 
au n" 26, a pour aire exacte 

F(a: + ft)— Fx = /tF' + iCT" + i/i'F"' +Lh'V"[x,,), 

La valeur approximative est ; 

— A[5F'(x)-f'8F'{3:-+- h) — V-{x + 'Zk)\ 

L'erreur est donc 

24 

nulle dans le cas d'une parabole du second degré, et, par suite 
du troisième dogré. On peut déduire ee résultat aussi du n° 32, 
111, et de la théorie de la parabole. 

Si l'on prend pour F'(«) l'ordonnée moyenne de l'aire consi- 
dérée, il est clair, d'après le n° 39, que l'erreur sera de la forme 

comme pour les deux formules noit\ellcs et pour eellcs de 
Simpson. 
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Comparaison nus formules, 

4S. Comparaison au point de vw de l'exacUtude. Si l'on géné- 
ralise les résultats trouvés dans le § Xli, en traitant divers exem- 
ples, et ceux que nous avons obtenus dans les g§ XVI et XVIII, 
moyennant les hypothèses du n" 36, on peut diviser les for- 
mules étudiées en trois groupes. 

Le premier comprend les deux formules de Simpson et celle 
de Weddie, qui sont exactes, les deux premières, pour les para- 
boles du S""" ou du 5"" degré, la dernière pour celles du 2""% 
du 3"% du 4™^ et du 3"°° degré. L'erreur pour les deux pre- 
mières est de la forme AhW(x.i), pour la dernière de la forme 
Bh'^F'"'(x/t). Pour arriver à cette estimation de l'erreur, on sup- 
pose, pour les premières, que F'(^)) F"(^)t ^"'(^^)' ^"(j^)' FX^) 
sont continues et que ¥Xx) varie peu, dans un intervalle 2A pour 
la première, dans un intervalle 3A pour la seconde. Pour la for- 
mule de Weddie, on admet que F'(sc) et ses six premières déri- 
vées sont continues dans un intervalle 6A et que F"" (a;) varie 
peu dans le même intervalle (*). A cause de ces hypothèses plus 
étendues et plus nombreuses, il est probable que la formule de 
Weddie ne peut pas être mise au-dessus de celles de Simpson. 
D'ailleurs, quand on compare les formules ici en question, d'après 
les limites supérieures de l'erreur commise, limites trouvées 
§g IV, X, XI, on reconnaît que la première formule de Simpson 
l'emporte sur les autres. On peut donc, en réunissant les divers 
éléments d'appréciation, conclure que les trois formules du pre- 
mier groupe ont à peu près la même valeur, au point de vue de 
l'exactitude. 



(*j (lue tonnulo obtenue par la méthode de Ç6les et esacta poiifutie parabole d'ùrdre n, 
conduit k une erreur de h forme A/i'-t-'F"!-* [x/x), pour une courbe ij = f'ix), moyennant 
(n + 1) conditions de coutiuiiité, dans un intervalle iili, et pourvu que F»+'(ji) varie peu 
dans cet intervalle. Le nombre des conditions pour que la forJuule soil très exacte croit 
donc très rapidement atec le degré d'exaoUtude. 
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Le second groupe contient lus formules de M. Catalan et de 
M. Dupain et la formule nouvelle. Elles sont exactes aussi pour 
les paraboles du 2'"' et du 3"" degré; l'erreur est de la forme 
Aft'F'(iK^); mais pour arriver à ce résultai, nous avons supposé 
que F'(x) et ses quatre premières dérivées sont continues dans 
un intervalle nh, où d'ailleurs P'(x) varie peu. A. cause de la 
grandeur de cet intervalle, en général, ces trois formules sont 
probablement moins exactes que les précédentes. 

Le troisième groupe comprend les autres formules. Celle de 
M. Parmentier donne nn résultat très voisin de celui que fournit 
la formule de M. Dupain; l'erreur ne dépend de h^ qu'à cause 
des éléments extrêmes, !l est donc extrêmement probable qu'elle 
vaut mieux, en général, que celle do Poneelet, et à fortiori que 
la formule des trapèzes (*). 

43. Comparaison au point de vue de la simplicité des calculs. 
Au point de vue de la simplicité des calculs, les neuf formules 

considérées forment encore trois groupes, presque identiques aux 
précédents. 

Le premier comprend la formule des trapèzes, qui est la plus 
facile à retenir, puis les formules de Poneelet et de M, Parmen- 
tier qui contiennent moins d'ordonnées, enfin la formule de 
M. Dupain, qui est la moins simple des quatre. 

La formule de M. Catalan et ia formule nouvelle forment le 
second groupe : elles sont un peu plus compliquées, à cause de 
la quantité d' qui y entre. 

Enfin, les formules dont l'application entraine les calculs les 
plus pénibles sont celles de Simpson et de Weddie, qui consti- 
tuent le troisième groupe. 

44. Conclusion. Si l'on n'a pas besoin d'une grande approxi- 
mation, le mieux, d'après ce qui précède, est de recourir à la 



(■) On irouve dans la Sialique graphique de Culmann, traduction l'raiiçaïse, Paris. 
Dunod.i880,pp. -103-103, une discusaïou sur la valeur comparée des formules de Simpson 
et de M, Pamientier au il essaie de prouver que la dci'nitre est supârieuro â la première. 
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formule des trapèzes, ou, si le nombre des ordonnées est impair, 
à celle de M. Parmentier, qui est plus exacte et exige le calcul 
de beaucoup moins d'ordonnées. 

Si Ton veut une exactitude plus grande, et si le nombre des 
ordonnées esl pair, on prendra la formule de M. Catalan. Quand 
le nombre des ordonnées esl impair, on se servira plutôt de la 
formule de M. Dupain , qui emploie moins d'ordonnées , ou , 
si l'exactitude esl requise avant tout, de la première formule de 
Simpson. 

Les autres formules peuvent être employées dans des cas spé- 
ciaux, par exemple, la seconde formule de Simpson, dans le cas 
où la base de l'aire à mesurer est divisée en 3« parties égales. 



H1STOBI1511E. 

4S. Formules de Simpson et de Weddle. Thomas Simpson (*) 
a publié sa formule dans ses Mathematical Dissertations, 1743, 
p. 109 ("). Au fond, il cherche l'aire des trapèzes curvilignes 
analogues à aABCc(Fig. 10), en substituant à la courbe ABC une 
parabole passant par les trois points A, B, C et dont l'axe est per- 
pendiculaire à ac. Cette démonstration a passé dans la plupart 
des traités de calcul intégral. — Poncclet, dans son Introduction 
à la Mécanique industrielle {n" 180, pp. 187-190 de la 2"° édi- 
tion, Metz, 1839; S-"» édition, Paris, 1870, pp. 194-197), en a 
donné une démonstration élémentaire comme il suit : Pour ti'ou- 



(*) Mé en 1710, mort en 4761; ne pas le coufondre avec son contemporain, Robert 
Sirason, né en 188T, mon en 1161 aussi. 

i") Hous empcunlons ce renseigoement a une note de la page 67, de la télÈbre édition 
de la Résistance des Corps solides de Havier, donnée par M. B. de Sainl-Vanant, Paris, 
Dunod, 1864. T. Simpson, dans son ouvrage inliwlé; The ùaea-hie and Application af 
Fluxions, LondonnSO, p. 202, lait allusion à sa formule, mais aeîa donne pas. Probable- 
ment, les Dissertations de Simpson contiennent aussi sa seconde formule, que nous avons 
emprunlée à un ouvrage de Rankine lUseful Riites and Tables. Foufth Edition. Londou. 
CGriffinandC-, 1813, p. 6S). 
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ver Taire nABCc, il divise oc en trois parties égales au, uv, vc, 
et mène les ordonnées de la courbe mU, uV, ainsi que les cordes 
AU, UV, VC, dont la seconde UV rencontre bB en W. L'aire 
curviligne nABCc est un peu plus grande que celle du polygone 
aAUVCc qui a pour mesure j h (j/^ -l- 4('W'-+- y^') ; en rempla- 
çant 46W par la quantité un peu plus grande iy^t ce qui revient 
à ajouter au polygone aAUVCc, quatre fois le petit triangle UBV, 
on se rapprochera vraisemblablement de l'aire exacte de aABCc. 
— MM. Rouehé et de Coraberousse, dans leur Traité de Géomé- 
trie (V" partie, paris, Gaulhier-Villars, 1879, pp. 294-296) ont 
simplifié cette démonstration en remarquant que si AC coupe ôB 
en b^, l'aire du trapèze aACc est égale à ^ h(y^ -t- i-bb^ -h i/j) ; 
remplaçant 466i par itj^, ce qui revient à ajouter [au trapèze les 
j du triangle ABC, on se rapprochera, vraisemblablement encore, 
de l'aire cherchée, — Il est clair qu'aucune de ces trois méthodes 
ne fournit une vraie démonstration de la formule , puisque l'er- 
reur maxima commise n'est pas évaluée (*). 




1") Dans les nouvelles Annales de Mulhèinatiqiies , [" série, t. XV, pp. 291-293, un 
Anonyme metlla tonniile de Simpson sous la forme s = ( 4- ^- (r — m'] at dît que Saygej 
en a donné une démonslration simple dans sa Géométrie élémemaire. Nous n'avons pn 
consulter cel ouvrage. 
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